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Abstract 

We re-examine excitation of a set of disk oscillations in a deformed disk by a 
resonant process. We assume that the disk is deformed from an axisymmetric steady 
state by an oscillatory deformation with frequency ud and azimuthal wavenumber 
TTT-D- Then, we consider two normal mode oscillations with a set of frequencies and 
azimuthal wavenumber being {ui, mi) and {002, and satisfying the resonant con- 
ditions {ui + + i^D = and mi + m2 + mo = 0). These oscillations are resonantly 
excited if E1E2 > 0, where Ei and E2 are wave energies of the above two oscillations, 
when the deformation is maintained by external forces or has a large amplitude com- 
pared with the oscillations. This instability condition is rather general as long as 
unperturbed density and pressure vanish on the surface of the system. Possibility of 
application to superhump and negative superhump in superoutburst state of dwarf 
novae are briefly discussed. 

Key words: accretion, accrection disks — disk deformation — instability — 
oscillations — resonance — superhump 



1. Introduction 

Kato et al. (2011) derived conditions of resonant excitation of disk oscillations in de- 
formed disks (see also Ferreira & Ogilvie 2008 and Oktariani et al. 2008 for numerical studies 
and Kato 2011 for additional analytical studies). The initial purpose of having examined this 
wave excitation process was to suggest a possible excitation mechanism of high-frequency quasi- 
periodic oscillations (QPOs) observed in black-hole low-mass X-ray binaries (BH LMXBs) (see 
Kato 2004 and subsequent papers, for example, Kato 2008, 2012). Not only our model but 
also many other models of HF QPOs in BH LMXBs have been proposed, but there is still no 
common consensus on the origin of the QPOs. 

To clarify whether the wave- wave resonant excitation process by Kato et al. (2011) is 
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really one of realistic processes of excitation of disk oscillations in deformed disks, it will be 
important to examine whether the mechanism can describe observations in systems where disks 
are definitely deformed and oscillations are observed. A typical example of such objects is dwarf 
novae. The gravitational fields of the disks of dwarf novae are tidally deformed by a secondary 
star, and if certain conditions are satisfied the disks become unstable (tidal instability) and 
oscillatory phenomena (superhumps) appears. Detailed observational studies on the oscillatory 
phenomena in dwarf novae have been performed by some groups, for example, by VSNET group 
by T. Kato (a recent work by his group is, for example, Ohshima et al. 2012). 

The tidal instability and superhump in dwarf navae are understood by the tidal instabil- 
ity (Lubow 1991) and precession of one-armed p-mode oscillation (Osaki 1985). The wave- wave 
resonant excitation process can also describe them (Kato 2013). In this contexts, we feel that 
the formulation by Kato et al. (2011) (hereafter paper I) should be improved, since i) not all 
possible resonant cases are considered there, ii) in some parts, nonlinear treatment of complex 
quantities was somewhat vague, and iii) interpretation of the instability was not enough. The 
purpose of this paper is thus to refine the analyses in paper I and to make the meaning of the 
instability more clearer. As an application, excitation of corrugation waves (warps) in dwarf 
novae disks, in addition to excitation of one-armed p-mode oscillations, is briefly mentioned. 

2. Outline and Assumptions in This Paper 



On an axisymmetric steady disk a deformation is superposed. The deformation is as- 
sumed to be time-periodic with frequency wd and with azimuthal wavenumber m^. In tidally 
deformed disks or in tilted disks, such deformation will be expected. Even in other situations, 
long-living time-periodic deformations will be present on disks by a long-scale instability of 
disks or by interactions with surrounding media. Our basic assumption here is that such a 
time-periodic deformation is maintained on disks over a timescale longer than the growing or 
damping timescale of oscillatory perturbations on disks. The displacement vector, ^Y){r,t), 
associated with the deformation over the steady state is denoted by 



^D(r)exp(iwDi) 



3ft 



(1) 



where 3ft denotes the real part, and (p is the azimuthal coordinate of the cylindrical coordinates 
(r, ip, z) whose center is at the disk center and the 2-axis is the rotating axis of the disk. 

In addition to such a deformation, a set of two normal modes of oscillations are super- 
posed on the axisymmetric steady disks. The set of eigen-frequency and azimuthal wavenumber 
of these oscillations are denoted by (wi, mi) and {002, 1712). The displacement vectors, ^i{r,t), 
associated with these oscillations are expressed as 



^i(r,t) = 3ft 



li(r)exp(iu;D^) 



3ft 



1,2). 



(2) 



We now assume that the following resonant conditions among the two oscillations and 
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the deformation are present: 

Ui + UJ2 + ujd = Au and mi + m2 + mD = 0, (3) 

where mi's (i = 1,2, and D) are integers. In order to include in our formulation the cases where 
three frequencies are slightly deviated from the exact resonant condition, Au is introduced in 
the first relation of equation ([3]), where |Aci;| is assumed to be much smaller than the absolute 
values of ui and uj2\j In the final analyses of instability, however, we consider only the case 
of Au = 0. In order to represent the resonant conditions by a simple form without separately 
considering such cases as uj2 = Ui + ud and uj2 = — cui +^05 we have adopted such resonant 
forms as equations ([3]), where w's and m's can be taken to be positive or negative as long as 
the conditions are satisfied. 

In the non-linear stages, the Ui- and a;2-oscillations satisfying conditions ([3]) resonantly 
interact with each other through the disk deformation specified by cjdS Our main concern in 
this paper is to examine how the amplitudes of the Ui- and a;2-oscillations are changed by the 
quasi-nonlinear interactions with the disk deformation. Quasi-nonlinear behavoirs of the Ui- 
and ci;2-oscillations themselves are not considered, assuming that their amplitudes are smaller 
than that of the disk deformation. It is further assumed that the amplitudes of the oscillations 
are so small that the amplitude of the disk deformation is unchanged by the resonant couplings 
with the oscillations. 

In order to understand the results concerning the instability condition, however, it is 
instructive to consider the cases where the amplitude of the disk deformation is not kept 
unchanged, but changes as results of resonant interactions with the oscillations. This is the case 
of three wave resonant interactions where no particular position is given to the disk deformation. 
This case is mentioned in section 5. Application to excitation of superhumps in dwarf novae is 
briefly mentioned in section 6. 

3. Formulation of Coupling Processes 

When we consider high-frequency QPOs in neutron-star and black-hole low-mass X-ray 
binaries, the effects of general relativity are important for quantitative studies. The effects 
of general relativity, however, is not essential in understanding the essence of the instability 
mechanism. Hence, in this paper, for simplicity, we formulate the problem in the Newtonian 
frame. We adopt a Lagrangian formulation introduced by Lynden-Bell and Ostriker (1987). 

^ The cases where wd = can be included in our analyses. 

^ For oscillations to resonantly interact, additional conditions concerning wave forms of the oscillations in the 
vertical and radial directions are necessary. These conditions are not mentioned here. If these consitions are 
not satisfied, the value of the coupling term, W or W"^ , given by equation (psi) or ([25)1 vanishes. 
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^°^-<5f-V^--Vp), (4) 



3.1. Hydrodynamical Equations for Linear Oscillations and Disk Deformation 

The unperturbed disk with no deformation is steady and axisymmetric. In the 
Lagrangian formulation, hydrodynamical perturbations, including a disk deformation, super- 
posed over the unperturbed disks can be described by (Lynden-Bell & Ostriker 1967) 

Dt^ "V ' p 

where ${r,t) is a displacement vector associated with the perturbations, and Dq/ Dt is the time 
derivative along an unperturbed flow, UQ{r), and is related to the Eulerian time derivative, 
d/dt, by 

In equation (jlj) 5{X) represents the Lagrangian variation of X, and if) is the gravitational 
potential. Other notations in equation (jlj) have their usual meanings. 

In the cases of small amplitude, adiabatic and inviscid perturbations equation (j4]) is 
written as (Lynden-Bell & Ostriker 1967) 

Po0 + 2poK-V)^ + L(O = O, (6) 

where L{$,) is a linear Hermitian operator with respect to ^. Equation ^ is valid even when 
the perturbations are self-gravitating. Hereafter, however, we consider only the cases of non- 
self gravitating perturbations. 

The disks where the behavoirs of oscillations are examined are not a steady equilibrium 
one, but in a deformed state as mentioned above. If the disk deformation is maintained by a 
tidal force, for example, the deformation is a forced oscillation and the displacement vector, 
^o, associated with the deformation is governed, in the linear approximations, by 

+ 2po(^o ■ V)^ + = -poV^D, (7) 

where ipD is the Eulerian variation of the gravitational potential, ip, due to the external tidal 
force, and is expressed as 



4jT){r)exp{iuY)t) 



(8) 



where cjd is the tidal frequency associated with azimuthal wavenumber m^, and cjd and 
are related by ud = mD^^orb, being the angular velocity of rotation of the secondary star. 

If no external force acts on the disk, the displacement vector, associated with the 
deformation is governed by an equation having the same form as equation ([6]) in the linear 
approximation. 

3.2. Quasi- Nonlinear Coupling among Oscillations and Disk Deformation 

Now, we consider two normal modes of oscillations, and with {ui, mi) and {u2, 
7712). In the linear stage, the perturbation, ^, imposed on a deformed disk is simply the sum of 
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these two oscillations: 



Ai$,^ex^[i{ujit - rriiLp)] + A2^2^xp[i{uJ2t - 1712^)] 



(9) 



where Ai and A2 are amplitudes and are arbitrary constants at the linear stage. In the quasi- 
nonlinear stage, the oscillations are modified through quasi-nonlinear couplings with disk de- 
formation, ^j-,, so that they satisfy a quasi-nonlinear wave equation. In the case where the 
deformation is internally maintained (i.e., no external force), the non-linear wave equation is 



Po 



2poK-V)^ + L(0=PoC(C,^d) 



(10) 



where C is the quasi-nonlinear coupling terms and is given by (Kato 2004, 2008) 



+ 



d 



Q3 



dridrjdr^ 



'0-2t^ POT^ TT^ 



dvj dri J 



d_ 

dri 



(ri 



l)Po 
l)Po 



f^div^. + f^div^ 



dvj drk 



d_ 

dri 



Ti{Ti - l)podiv| ■ div^D 



(11) 



where the subscript i attached to C and r denotes the i-component of vectors C and r in the 
Cartesian coordinates. Furthermore, Fi is the barotropic index specifying the linear part of the 
relation between the Lagrangian variations of pressure, 5p, and density, 5p, i.e., 5p/po = ^i^p/ Po- 
In the cases where the disk deformation is due to an external force, additional terms 
should be added to equation (ITTl) . In the case of tidal deformation, for example, the Eulerian 
variation of the gravitational potential, ipD, is associated with ^y)- Because of this, in con- 
sidering the Lagrangian variation of Vip [see SVip in equation (j4])], a quasi-nonlinear term of 
^jdVipB/drj, which has been neglected in deriving equation ( ITT]) , should be taken into account, 
since we have 



d 1 d'^ 

VV'D + 07^V(^o + ^d) + i^^jCkT^rnr^i^o + ^d) + 



' drjdrk 



(12) 



where iJjq is the unperturbed gravitational potential in the unperturbed axisymmetric disk. 
Corresponding to this, the coupling term C is changed to C^, which is 

C^it^B) = C(^,|d) -e.T^V^D. (13) 



drj 



Here and hereafter, the superscript T is attached to C, when the disk deformation is due to 
tidal force, and the quasi-nonlinear wave equation is 



Po 



(14) 



The disk oscillations, ^(r,t), resulting from the quasi-nonlinear coupling through disk 
deformation, ^dI^;^); '^i^ be written generally in the formjf] 



The expression for ^{r,t) in Kato et al. (2011) was somewhat vague. 
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2 2 
i=l i 07^1,2 

+oscillating terms with other frequencies. 



(15) 



The original two oscillations, and resonantly interact through the disk deformation. 
Hence, their amplitudes secularly change with time, which is taken into account in equation 
(fT5l) by taking the amplitudes, A^s, to be slowly varying functions of time. The spatial forms 
of the terms whose time- dependence is exp{iuit) but whose spatial dependence is different from 
are expressed by a sum of a series of eigen-functions, (« 7^ 1 and 2), assuming that they 
make a complete set. The terms whose time-dependences are different from exp[iuit) are not 
written down explicitly in equation (fT5l) . since these terms disappear by taking long-term time 
average when we are interested in phenomena with frequencies of cui (i = 1 and 2). 

In order to derive equations describing the time evolution of Ai(t), we substitute equation 
(fT5i) into the left-hand side of equation (fTOl) or ( |T4l) . Then, considering that ^j's and ^^.'s are 
displacement vectors associated with eigen-functions of the linear wave equation ([6]), we find 
that the left-hand side of equation ( |TOl) or (fT4l) becomes the real part of 



i=l 



dt 



[uo ■ V) 



liexp(iwit) 

- 2i{ljJa — COi)iUo ■ V) 



(16) 



where d?Ai/dt^ has been neglected, since Ai{t){i = 1 and 2) are slowly varying functions of 



time. Now, the real part of equation (fT6 
being multiplied by £i[= 3? $iexp{iuit)] 



is integrated over the whole volume of disks after 
Then, the term resulting from the second term of 
equation ( fT6i) vanishes [see an orthogonal relation given by equation (8) of paper I] and the 
results become 

.dAi 



^ i- 



dt 



(17) 



where (X) denotes the volume integration of X{r), and the asterisk * denotes the complex 
conjugate. This can be further reduced to 

* .^-^'-^ (18) 

ui dt ^ ^ 

by using the wave energy of the oscillation of ^i, Ei, which is a real quantity defined by (e.g., 
Kato 2001) 

1 



El 



wi(Po^I^i) -«(po^I(*to- V)^i) 



(19) 



1, 



1, 



The formula 

5R(A)SR(B) + AB*] ^ ^^AB + A*B] 

is used, where A and B arc complex variables and B* is the complex conjugate of B. 
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It is of importance to note that the wave energy can be expressed in the case of the oscillations 
in geometrically thin disks as (e.g., see Kato 2001) 

E, = ^(^{co, - rmn)poiilM,r + ilAz)), (20) 

where Q{r) is the angular velocity of disk rotation. 

Now, the real part of equation ffTOj) is multiplied by $,i{r) and integrated over the whole 
volume. Then, considering the results of the above paragraph, we have 



Ui at 2 



A2{t)AJi, ■ C(^2, ^d) )expitALot) 



(21) 



where C(^2)^d) is replaced to C^(^27^d)) when tidal deformation is considered. In deriving 
the right-hand side of equation fl^Tl) . the time periodic terms with high frequencies (i.e., non- 
resonant terms) have been neglected by time average being taken. 

Similarly, we multiply ^2(^) ^he real part of equation ffTOl) and integrate over the 
whole volume to lead to 

"^'^^ = 1^ [A(t)^D(i-C(^i,CD))exp(^Aa;t)]. (22) 

As mentioned above, in the case of tidal deformation C(^]^,^d) is replaced by C^(^i,^d). 

4. Growth of Resonant Oscillations 

Since the amplitude is assumed to be constant, the time evolutions of Ai and A2 
are determined by solving the simultaneous equations ( 1211) and ( l22l) . What are governed by 
equations fl2T|) and fl22l) are the imaginary part of Ai and A2, i.e., A^ and A2i. Their real 
parts are not related to the resonance, and we can neglect them in considering the exponential 
growth (or damping) of An and A2i. Hereafter, we restrict our attention only to the case of 
Au = 0. Then, we have from equations (12T!) and (l22l) 

^^'"^'^ Ia2MAbW), (23) 



oji dt 2 

2^2 A2,i _ 1 

UJ2 dt 2 



A,^^{A^W), (24) 



W^i,- C(i,^D) ) =(^2 ■ C(^i,Cd) )• (25) 



where 



The last equality showing commutability of and ^2 is important for obtaining an expression 
for instability. This commutative relation is derived by using an expression for C given by 
equation (ITTi) and performing the volume integration by part, neglecting surface integrals under 
the assumption that the unperturbed density and pressure vanish on the surface of the system. 
In the case where the deformation is due to the tidal force, W^s in equations (1231) and flM|) are 
replaced by which is 



(26) 



The commutability in equation (!26|) is obtained by using given by equation (fT3|) . 



Eliminating A2,i from equations (l23i) and ( l24l) . we have 




(27) 



The same equation is obtained for by ehminating Ai^i instead of A2;i. 



Equation ( 1271) shows that if 




(28) 



the oscillations grow with the growth rate given by 



(29) 



In the case where the disk deformation is due to tidal force, W in equation (129|1 is changed to 



At a glance, the sign of the instability condition given by equation fl28|) is opposite to 
that in paper I. This is due to a formal difference of expression for the resonant condition, 
and there is no difference in contents. In paper I we have examined two resonant cases of 
ul =ujI + uj\^ and ul = — (the superscript I is attached to w's in order to distinct these w's 
from those in the resonant condition of + u;2 + wd = in the present paper), and derived the 
instability condition: EiE2/u}\ul < 0. The resonant case of ul = uj\+uj\^ in paper I corresponds 
to the case where 002 in the present paper is formally written as —uj^. If this change is adopted 
the present instability condition becomes EiE2/u)iUj\ < 0, and is the same as that in paper I. 
Furthermore, the case of u\ = — in paper I corresponds to the case where Ui in the present 
paper is written as —uj\ and we have EiE2/uj\uj2 < 0. 

In principle, the cases of = —{uj\ + w^) and ul = — (w{ — w^) are possible, but not 
examined in paper I. If we consider these cases by the same way as in paper I, we obtain 
EiE2/u}\ul > 0, different from the cases of ul = u\+ujI) and ul = u\ — This instability 
condition is also consistent with the condition, E1E2/ UJ1UJ2 > 0, of the present paper. This is 
because in these cases ui = u\ and 002 = oul (the sign of ub does not appear in the instability 
criterion). In summary, in all resonant cases the instability condition is summaried by equation 
(j28l) . if the resonant condition is written as Ui + UJ2 + ujb = ^■ 

5. Three Wave Interactions and Instability 

In the case where \ud\ is small, the resonance occurs between two oscillations with 
opposite signs of w's, i.e., ujiuj2 < 0. In this case the instability condition (l28l) is written as 



E1E2 < 0. 



(30) 
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ation is 

i 



This might suggest that the instabihty is, at the lowest order of approximation, a resuh of energy 
exchange between two oscillations through a catalytic action of disk deformation. Positive 
energy flows from a negative-energy oscillation to a positive-energy oscillation, leading to growth 
of both oscillations. If cji > and uj2 <0 (and also mi > and m2 < 0), i.e., if both oscillations 
are prograde with opposite signs of frequency, we introduce U2 and m2 defined, respectively, by 
—U2 and — 777,2 so that all w's and m's become positive. Then, the instability condition ( l30l) is 
reduced to 

{ui — mifl) {0J2 — 1712^1) < 0, (31) 

where ( X ) is a weighted mean of X in the region where oscillations prevails [see expression 
(I20I) for wave energy]. This equation shows that for the instability to occur, one oscil 
inside of its corotation radius, while the other must be outside of its corotation radius 

The cause of instability is, however, not always interpreted in terms of energy exchange 
between ui- and a;2-oscillations alone. Energy exchange between deformation and oscillation 
is of importance in the case where \ujy)\ is large compared with and \uj2\^ In this case the 
resonance {ui + uj2 + u}d = 0) occurs for oscillations with the same signs of ui and U2, and the 
instability condition ( 128|) is written as 

£^1^2 > 0, (32) 

or 

(ui — miQ) {uj2 — rn2^^) > 0, (33) 

distinct from the case where I^dI is small. Equation (l32i) clearly shows that in the present 
case of \ujy)\ being large, the instability cannot be interpreted as energy exchange between two 
oscillations. That is, both oscillations grow by getting energy from the deformation (in the 
case of i?i > and £'2 > 0), or by lossing energy to the deformation (in the case of -Ei < and 
£"2 < 0). In other words, energy exchange between oscillations and deformation is esential for 
growth of the oscillations. 

Unless the disk deformation is maintained like the case of tidal deformation, or the 
deformation has a sufficiently large amplitude, the effects of amphtude change of the disk de- 
formation on oscillations will be non-negligible. To examine this, let us consider the case where 
the deformation is not maintained, but simply follows the same hydrodymnamical equation as 
those of the oscillations. Then, by similar processes with those used to derive equation (12T!) 
and fl22l). we have 



(34) 



The corotation radius of an oscillation with uj and m is defined as the radius where w = mO. Inside the 
radius uj — niQ < 0, while outside the radius uj — mQ > 0. 

The tidal instability in dwarf novae corresponds to this case. In paper I this case was not explicitly 
considered. 
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where Ed is the wave energy of disk deformation. In this case we have the commutative relation: 
'^D ■ C{ki,L)) =/i ■ C(|2,Id)) =fe ■ C(i,lD)\ = W. (35) 



Hence, if Au = 0, we have the following set of equations [see equations fl23|l and ([2 
2Ei dAi i 1 



Ui dt 2 



A2,iAz),i3fiW^, (36) 



-U,,ADm. (37) 



UJ2 dt 2 

'""^^^^^^ -\A,,A2m- (38) 



dt 2 

In order to examine how the amplitude, Ao^iit), changes with time, we take the time 
derivative of equation ( l38l) and substitute equations ( l36i) and ( 137|) into the resulting equation. 
Then we have 



(Siy)2Az),i. (39) 



(it^ 16Vu;2u;d 

This equation shows that in the case of {Ei/uji){E2/uj2) > 0, Ao,i damps with time unless 
[En /ujd){E2/u)2) > 0. However, {Ei/uji){E2/uj2) > and {E£,/ojd){E2/uj2) > cannot be realized 
usually at the same time in resonance. This comes from the following situation. The sign of 
Ed/ujd is equal to the sign of (uD — rrio^), where ( X ) is a some spatial average of X [remember 
the definintion of wave energy, i.e., see equation (120|) ]. Hence, when the resonant conditions are 
satisfied, we have 

Sign {ojo — mjjVL) = —Sign [{ui — rriiQ) + {uj2 — 1712^1)]. (40) 

Roughly speaking, the signs of {ui —miQ) and {uj2 — m2^l) are equal, respectively, to the signs 
of Ei/uji and E2/UJ2, since for resonance to occur the propagation regions of all the oscillations 
must overlapp. This means 

Sign — < 0, if — > and — > 0, (41) 

COd UJ2 

Sign — > 0, if — < and — < 0. (42) 

OJD Wi UJ2 

In both cases we have {Eij/ujo){E2/uj2) < 0. In other words, when Ui- and aj2-oscillations grow 
by (-El /oJi) (-£^2/1^2) > 0, WD-oscillation damps. 

This situation can be also demonstrated clearly from the following considerations. 
Equations fl36l) - fl38|l are a set of simultaneous non-linear differential equations with respect to 
^2,i, and Ad^i- By multiplying i and ^42,1, respectively, to equation fl36|) and fl37j) . and 
taking the difference of the resulting two equations, we have 

—Al^-—Al^= const. (43) 

Ul ' UJ2 ' 

By the similar procedures, we also have 
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ElAl.-^Al,= const. (44) 
^Al,-^Al,= const., (45) 

U2 ' CiJd 



where the right-hand sides of equations ( 1431) - fl45|) are time-independent constants. Equation 
( l43l) shows that when {Ei/ui){E2/oJ2) > 0, Ai^i and A2,i can increase with time. Equation (jS]), 
however, demonstrates that in that case A^^i decreases since the sign of E2/U2 and E^/ud is 
opposite. 

More generally speaking, let us consider resonant interactions among three oscillations. 
Then, two of them can grow by consumption of the other oscillation, but not all oscillations 
can grow. 

6. Applications to Superhumps in Dwarf Novae 

Typical objects with deformed disks which are maintained are dwarf novae, where the 
gravitational potential in disks is tidally deformed by a secondary star. Whitehurst (1988) found 
that under a certain condition the disks become unstable to an eccentric deformation having 
slow precession. The instability is known to be due to a tidal instability of the so-called 3 : 1 
resonance (Lubow 1991). The precession comes from excitation of one-armed oscillation and it 
brings about superhump phenomena (Osaki 1985). (See Osaki 1996 for review of superoutburst 
and superhump of dwarf novae.) Lubow (1992) further shows that a tilt is also excited at the 
3 : 1 resonance by a mode-coupling process similar to that in the case of excitation of the 
eccentric deformation. The present wave-wave resonant instability process can also describe 
the instability and excitations of both one-armed oscillation and tilt at the 3 : 1 resonance. 
Although the excitation of one-armed oscillation was already discussed in Kato (2013), we 
briefly summarize here the essence of excitation of both modes. 

6.1. Tidal Instability and Excitation of Low- Frequency One-Armed p-Mode Oscillation 

In tidally deformed disks, mo and ud are related by ud = mofiorb; where Qorh is the 
orbital frequency of the secondary star and niD is an integer. As the wi-oscillation, we take 
one-armed (mi = 1) low-frequency p-mode oscillation with a positive frequency. The dispersion 
relatior0 shows that the frequency, cui, must be smaller than f2 — k, i.e., wi < fi — k, where Q{r) is 
the angular velocity of disk rotation and K(r) is the epicyclic frequency in the equatorial plane. 



Let us consider vertically isothermal disks. Then, the dispersion relation for perturbations which are local 
in the radial direction is given by 

where n± is the vertical epicyclic frequency, n an integer representing the node number of oscillations in the 
vertical direction, k the radial wavenumber, and Cs the acoustic speed in disks (see, e.g., Kato 2001 or Kato 
et al. 2008). It is noticed that ^± > il> k in geometrically thin disks. 
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Fig. 1. Schematic diagram showing frequencies and propagation regions of the wi- and a;2-oscillations in 
the case where the wi -oscillation is one-armed low- frequency p-mode oscillation. The scales of coordinates 
are arbitrary, and are not linear. The one-armed p-mode oscillation is trapped between Tc and r^. The 
inside is the evanescent region. The a;2 -oscillation can propagate inside r^. See figure 2 in Kato (2013), but 
notice that the signs of uj2 and m2 in this paper are adopted to be opposite to those in paper I. Because 
of this change, all the curves related to W2-oscillations in this paper are line-symmetric with respect to the 
horizontal line of w = to those of figure 2 in Kato (2013). 

Since VL{r) — K(r) is positive and increases outwards in tidally deformed disks, the oscillation 
with frequency Ui is trapped between the radius, Tc, where Wi = (fi — k)c (the subscript c denotes 
the value at r^) and the radius, rt(> r^), where the disk is truncated by the tidal effects (see 
figure 1). This oscillation has Ei/ui < 0, since ui — miVL = wi — < in the region where the 
oscillation dominantly exists. 

Next, let us consider a p-mode oscillation with U2 and m2. Here, UJ2 and m2 are taken 
to be positive, and m2 is found later to be 2, i.e., m2 = 2. The dispersion relation shows that 
one of the propagation region of the oscillation is specified by 072 — 'm2Vt < —k. Since 2VL — k 
is positive and decreases outwards, the propagation region of the oscillation is bound outside. 
The outer edge of the propagation region is taken to be (see figure 1) (see Kato 2013 for 
more discussion). This oscillation has E2/0J2 < 0, since U2 — m2^l < 0. 

In summary, we have (-Ei/cJi) (£'2/1^2) > 0. Furthermore, the propagation regions of the 
Ui- and a;2-oscillations are overlapped. Thus, the above set of oscillations are excited if the set 
of resonant conditions, ui + U2 + uj-d = and mi + m2 + niB = 0, are realized inside the disk. 
These resonant conditions are summarized as 



[1 + mr,)^ 



orb 



0, 



(46) 



where mi = 1 and m2 = — (1 + ^d) are used. For simplicity, if the difference between Q and k 
is neglected, the above equation leads to 



fit 



mD 



orb- 



(47) 
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Fig. 2. Schematic diagram showing frequencies and propagation regions of the ui- and a;2-oscillations 
in the case where the wi-oscillation is one-armed low-frequency corrugation mode. The scales of coordi- 
nates are arbitrary, and are not linear. The corrugation wave can propagates all region inside r^. The 
a;2-oscillation can also propagate inside r^. 

For m-Q = — the resonant condition is satisfied by taking so that fit '■ ^orb = 3 : 1 is 
reahzed. This can be really expected inside the disk. The one-armed oscillation excited is the 
cause of superhump (Osaki 1985). 

6.2. Excitation of Low- Frequency Corrugation Mode 

Let us show that an another set of oscillations is excited in the case where fit '■ fiorb = 
3 : 1 is realized. We consider the corrugation mode (i.e., the mode of tilt or warp with a 
slow frequency) as the wi-oscillation (for classification of oscillation modes in disks, see, for 
example, Kato 2001 and Kato et al. 2008). The mode is characterized by mi = 1 and ni = 1, 
where Ui is the node number in the vertical direction of the radial component of displacement 
vector, ^r{i^,t). The vertical component of the displacement, C,z{f',t), has no node in the vertical 
direction in the case of n = 1, and the disk plane at a given radius oscillates up and down in 
the vertical direction with the same phase. 

The dispersion relation shows that the propagation region of the low-frequency corru- 
gation mode is specified by — fi < — fi^, where fi^(r) is the vertical epicyclic frequency, and 
is larger than fi(r) in the tidally deformed disks, i.e., fi^ > fi, since the restoring force in the 
vertical direction increases by the gravitational force of the secondary star. Since fi — fi^ < 
and its absolute value increases outwards to the secondary star, the propagation region of the 
oscillation is inside the radius rt' specified by = (fi — fi_L)t', where the subscript t' denotes 
the value at rf. The frequency of the oscillation is negative, i.e., Ui < (see figure 2). This 
oscillation has Ei/ui < 0, since in the propagation region wi — fi < 0. 

^ It is noted that the formal expression for resonance in this paper, i.e., lui +uj2 +wd = 0, is different from 
that of Kato (2013), i.e., uj2 = ioi — wd- This causes differences of values of w's and m's in this paper from 
those in Kato (2013), but this is formal and there is no difference in real results. 
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As the (X'2-oscillation we adopt here the vertical p-mode oscillation with m2 = 2 and 
n2 = 1. So far, we did not mention about a resonant condition among ni, n2, and ud- Unless 
they satisfy a certain relation, however, the coupling term, C, vanishes and there is no resonant 
coupling. In this paper a plane-symmetric tidal force is implicitly considered, which corresponds 
to riD = 0. Hence, for resonant coupling to occur, the a;2-oscillation must have n2 = 1 in the 
vertically isothermal disksj^ The dispersion relation shows that one of the propagation region of 
the vertical p-mode oscillation with n2 = 1 is specified by U2 — 'm2^l < —Q±. Now, we take m2 = 2. 
Then, since 2Q — Qj_ is positive and decreases outward, the propagation region of oscillation 
with U2 is inside the radius rf specified by uj2 = (2fi — f2j_)t" > 0, where the subscript t" denotes 
the value at rt". This oscillation has E2/UJ2 < 0, since in the propagation region U2 — m2^l is 
negative. So far, the relation between rt' and rt" is not mentioned. We suppose that rt' ~ rt", 
but careful considerations are necessary to have a reliable relation. Here, we proceed without 
considering this problem in detail. 

Since {Ei/ ui){E2/ ^2) > 0, the above set of oscillations is excited if their propagation 
regions are spatially overlapped and the resonant conditions, ui + UJ2 + ujtj = and mi + 1712 + 
rriB = 0, are satisfied in the overlapped region. The resonant conditions are summarized as 

{n - fix)t' + (21^ - ^±)t" + mo^^orb = 0. (48) 
The radius determined by this relation depends little on rt' and we have 

fit" + JTlofiorb = 0. (49) 

This is satisfied for = ~3 at the radius of fit" • fiorb = 3 : 1. In summary, we have 
Wi ~ (fi — fi_L)t' < with mi = 1 and U2 ~ (2fi — fi_L)t" > with m2 = 2. The signs of Ui and 
mi are opposite, meaning that the oscillation is a retrograde wave. 

It is interesting that oscillations which slowly progrades (one-armed p-mode oscillation) 
and slowly retrogrades (corrugation wave) are both excited when the 3 : 1 resonance is realized. 
In dwarf novae, negative superhumps are often observed in addition to (positive) superhumps. 
The period of the negative superhumps is slightly shorter than the orbital period, different 
from those of (positive) superhumps. Furthermore, different from (positive) superhumps which 
are observed only in superoutburst, negative superhumps are observed in superoutburst (e.g., 
Ohshima et al. 2012), as well as in quiescent. The standard interpretation of the nagative 
superhumps is that the disks are tilted with retrograde precession and the gas stream from a 
secondary star hits the different part of the tilted disks (Wood et al. 2011). By analyzing in 
detail the observational data of V1504 Cygni by the Kepler, Osaki and Kato (2013) show that 
the light curve of the negative superhump in superoutburst is different from that in quiescent 
state, and suggest a model concerning the origin of the negative superhump in superoutburst. 

^ In the vertically isothermal disks, the 2-dependence of displacement vectors is expressed by use of the 
Hcrmite polynomials, and the orthogonality of oscillation modes are classified by the index of the Hcrmitc 
polynomials. 
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Our results that corrugation waves, which slowly retrograde, can be excited if the 3 : 1 resonance 
is realized suggest that the negative superhumps in superoutburst may be related to excitation 
of the corrugation waves. 



7. Discussion 

We have shown that in deformed disks where the deformation is maintained by some 
external force or it has a large amplitude compared with those of disk oscillations, a set of 
oscillations satisfying 



are excited on the disks, if they are eigen-mode oscillations in disks and satisfy resonant condi- 
tions, i.e., equations ([3]). This problem was already examined by paper I, but there were some 
points to be improved in their formulation and in interpretation of the results. This is the reason 
why we revisited to this problem. The original purpose of considering this excitation mecha- 
nism was to suggest a possible cause of high-frequency quasi-periodic oscillations (HF QPOs) 
observed in low- mass X-ray binaries (LMXBs). However, the objects which definitely have de- 
formed disks and where oscillatory phenomena are really observed are dwarf novae. Hence, as 
an example of applications of the resonant excitation process, we have demonstrated in section 
6 how the wave-wave resonant processes can describe the tidal instability and superhump in 
dwarf novae. 

The reason why the instability condition has such a simple form as equation (!50|) comes 
from the commutative relation of the coupling terms, i.e., equations ( l25l) and ( 126|) . These 
relations are quite general. Even in hydromagnetic systems these relations exist, which will be 
shown in a subsequent paper. 

It should be noticed that for uji- and a;2-oscillations to grow, the disk deformation must 
be externally maintained as in the case of tidal deformation, or has a sufficiently large amplitude 
compared with the oscillations. The HF QPOs in black-hole LMXBs appear only in the very 
high state (the steep power-law state). If the wave- wave resonant process considered here is the 
cause of HF QPOs, it suggest that the disks in the very high state are deformed from steady 
axisjTumetric state. Really, there are numerical simulations which show that in the very high 
state a thick torus is formed in the innermost part of the disk and a spiral pattern appears 
there (Machida & Matsumoto 2008). 

The auther thanks Y. Osaki for valuable discussions on negative superhump and A.T. 
Okazaki for discussions at the stage where paper I was written. 
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